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Abstract 

In three dimensional Fermi gases with spin imbalance, a competition exists between Cooper 
pairing with zero and with finite momentum. The latter gives rise to the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) superfluid phase, which only exists in a restricted area of the phase diagram as 
a function of chemical potential imbalance and interaction strength. Applying an optical potential 
along one direction enhances the FFLO region in this phase diagram. In this paper, we construct 
the phase diagram as a function of polarization and interaction strength in order to study the 
competition between the FFLO phase and the spin balanced BCS phase. This allows to take into 
account the region of phase separation, and provides a more direct connection with experiment. 
Subsequently, we investigate the effects of the wavelength and the depth of the optical potential, 
which is applied along one direction, on the FFLO state. It is shown that the FFLO state can exist 
up to a higher level of spin imbalance if the wavelength of the optical potential becomes smaller. 
Our results give rise to an interesting effect: the maximal polarization at which the FFLO state can 
exist, decreases when the interaction strength exceeds a certain critical value. This counterintuitive 
phenomenon is discussed and the connection to the optical potential is explained. 
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I. INTRODUCTION 



In the past decennium, the field of cold atoms has witnessed the experimental realization 
of various new quantum coherent phenomena in ultracold gases [1]. Amongst the various 
systems that are being studied, fermionic superfluids occupy a prominent role 2^4]. Con- 
trary to bosonic gases, fermionic gases must form pairs between fermions in different spin 
states in order to form an s-wave superfluid. A fundamental question that has attracted 
wide attention recently, is what happens when the ratio between these different spin-states 
is altered, thereby creating a spin-imbalanced or polarized superfluid. This effect is akin 
to applying a magnetic field to a superconductor, which leads to electrons aligning their 
spins. In superconductors, and in condensed matter systems in general, it is difficult to 
control various parameters such as the polarization or the interaction strength. In ultracold 
gases, however, one has an unprecedented degree of control over these parameters 
For instance, experimentalists can tune the interaction strength by making use of Feshbach 
resonances (for a review on this subject see Ref. 
from a Bardeen-Cooper-Schrieffer (BCS^ 
of bound molecules has been studied (9 



8|). Using this technique, the crossover 
superfluid to a Bose-Einstein condensate (BEC) 
121 ] . Furthermore, polarization can be created by 



first loading particles into one hyperfine state and subsequently applying a radio-frequency 
sweep to send a controlled number of particles to a different hyperfine state. This experi- 
mental liberty has led to the realization of many experiments on the imbalanced ultracold 
Fermi gas 13KL7I]. Today, the main phases of this system, as a function of temperature, 
interaction strength and polarization, have been observed experimentally 19jj . One of the 
most fundamental effects that occurs in a polarized Fermi gas is that above a critical im- 
balance, fermionic superfluidity ceases to exist and the system makes a transition into a 
normal Fermi gas. This is known as the Clogston-Chandrasekhar limit 18J. One important 
question that has emerged is whether there exist other pairing-mechanisms by which a gas 
of interacting fermions can accommodate polarized superfluidity. Recently, there has been 



an intensive theoretical 
or Sarma phase 



search for these exotic superfluid phases, such as the breached pair 



21] , and phase separation [22j. Another phase that has attracted wide 



attention is the Fulde-Ferrell-Larkin-Ovchinnikov phase (FFLO phase), that was proposed 



independently by Fulde and Ferrell 



23| and by Larkin and Ovchinnikov |24j in 1964. Their 



idea was that a fermionic system can accommodate spin imbalance, by forming pairs with 
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finite center-of-mass momentum. Up till now, this state has never been observed. A recent 
paper reports indirect experimental evidence for the FFLO state in one dimension (ID) 25], 
but in three dimensions (3D), the FFLO state has so far eluded experimental observation. 
This can be related to theoretical predictions, which state that the FFLO state only occupies 
a tiny sliver of the BCS-BEC crossover phase diagram It was therefore necessary 



to find a new method to stabilize the FFLO state in 3D. At the moment, new ideas are 
emerging to realize this goal. For instance, in two recent pap ers, it was suggested to stabi- 



lize the FFLO state by the use of a 3D cubic lattice |28|, |29|]. As an alternative approach, 
we proposed to stabilize the FFLO state in an imbalanced 3D Fermi gas by adding a ID 
optical potential to this system 30|. This turned out to enhance the stability region of the 



FFLO state by up to a factor six. In Ref. 



301 ]. the proof of principle for this concept was 



established using the path integral method , but questions remained about the exact effects 
of the properties of the optical potential on the FFLO state. In this paper, we investigate 
the role of the optical potential in detail. Our aim here is to present results that relate 
closely to the experiment. To this end, we construct and discuss the phase diagrams of an 
imbalanced Fermi gas in 3D, subjected to a ID optical potential, both at fixed chemical 
potentials and at fixed densities. Subsequently, we investigate and discuss the effects of the 
properties of the optical potential on the BCS superfluid state and on the FFLO state. In 
Sec. [Til we outline our strategy to construct the phase diagrams at fixed chemical potentials 
and at fixed densities. Subsequently, we discuss the competition between the BCS phase and 
the FFLO phase in these phase diagrams, which differ qualitatively from the corresponding 
phase diagrams for the imbalanced Fermi gas in 3D without the optical potential. In Sec. 
PU we investigate the effects of the parameters of the optical potential on the FFLO state. 
Furthermore, we explain the remarkable effect that the maximal polarization at which the 
FFLO state can exist decreases when the interaction strength is increased above a certain 
critical value. Finally in section IIVI we draw conclusions. 
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II. CONSTRUCTION OF THE PHASE DIAGRAMS 



A. Phase diagram at fixed chemical potentials 



To construct the phase diagram of an imbalanced 3D Fermi gas, subjected to a ID optical 
potential, at fixed chemical potentials, we start from the saddle-point free energy per unit 
volume of this system, which is given by jsoj]: 
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To derive expression ([T]), the ID optical potential was described by using a modified disper- 
sion relation for the fermionic particles: 
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where 5 is a function of the depth of the potential V and of the recoil energy E R , and Ql is 
the wave vector of the optical potential. Expression ([3]) consists of a free particle dispersion in 
ihe directions perpendicular to the optical potential, and a tight-binding periodic dispersion 



32| in the direction parallel to the optical potential. Because we want to describe the 



FFLO state, we only consider the BCS-side of the BCS-BEC crossover [33|, [34|. In this 
region, the scattering length is always negative, hence the following units are used ([Q): 
h = 2m = —a s = 1, where a s is the s-wave scattering length. Here it must be mentioned 
that the derivation leading to the free energy ([1]) neglected the density modulation of the 
Fermi gas due to the ID optical potential. This approximation leads to an underestimation 
of the interactions between the spin-up and spin-down fermions. 

The free energy ([1]) depends on two thermodynamic variables fi and (, and on two 
variational parameters A en Q. The variables \x and ( are the total chemical poten- 
tial jj, = (m + /ij /2, which fixes the total density, and the imbalance chemical potential 
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( = (fif — fii) /2, which fixes the polarization. A description in terms of \x and ( is equiva- 
lent to a description in terms of the respective chemical potentials of both spin species fj,^ 
and The variational parameters A and Q denote the binding energy and the center- 
of-mass momentum of the fermion pairs respectively. The parameter Q was introduced to 
include the FFLO state (which is characterized by fermion pairs with finite momentum) in 
our description. To find the ground state of the system, given fixed values of /i and (, the 
saddle-point free energy Q sp (/i, £; A, Q) has to be minimized with respect to the variational 
parameters A and Q. This defines the saddle-point equations 



which have to be satisfied simultaneously. One must be careful when extremizing the free 
energy, because only the minima correspond to a stable physical state. We avoided this 
difficulty by explicitly minimizing the free energy, instead of searching for roots of (jlj). 
Depending on the values of the parameters /i and (, three different ground states can emerge: 
1) a spin-balanced (S-b) superfluid, with A ^ and Q = , 2) the FFLO state with 
A ^ and Q ^ 0, and 3) the normal state with A = 0. Moreover, in the phase diagrams 
we will make a distinction between the partially polarized normal state (PPN) and the 
fully polarized normal state (FPN). In the latter, all particles are in the same spin-state, 
while in the former both spin-states are still present. Following this procedure, the phase 
diagram as a function of /i and ( can be constructed. Before doing so, the interval of 
the chemical potential fi that is relevant to describe the FFLO state must be estimated. 
As mentioned earlier, we only consider the BCS-side of the BCS-BEC crossover. This 
means that the variable (kpag)^ 1 , with kp the Fermi wave vector, which characterizes the 
interaction strength of a two-component Fermi gas in 3D, has to lie in the interval ]— oo, 0]. 
At unitarity, {kpa^ 1 = 0, which means that, in our description, the density has to go to 
infinity since we work in units where a s = —1. Even at {kfd s )~ l = 0, the chemical potential 
is still a fraction of the Fermi energy, and since the Fermi energy goes to infinity it follows 
that fi = +oo. In practice, it turns out that [i ~ 6 lies close enough to unitarity to describe 
the full stability region of the FFLO phase. The BCS-limit is already attained within good 
approximation at (^a^)" 1 = —2. In this limit, /i = Ep from which it follows that /i = 1/4. 
Figure [1] shows a comparison between the phase diagrams in the (fi, £)-plane of an imbal- 
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anced Fermi gas in 3D (figure [I] (A)), and an imbalanced Fermi gas in 3D, subjected to a ID 
optical potential (figure CD (B)). In both these phase diagrams the ground state of the system 




FIG. 1: Comparison between the phase diagram in the (/i, £)-plane (^ and £ are in units Eq = 2 ma 1 ) 
of (A) an imbalanced Fermi gas in 3D, and (B) an imbalanced Fermi gas in 3D, subjected to a 
ID optical potential. Both phase diagrams are constructed based on a finite set of points in the 
(//, £)-plane. For each point, the ground state of the system is determined. There are four possible 
ground states: a spin-balanced (S-b) superfluid (red circles), the FFLO phase (green triangles), 
the partially polarized normal (PPN) phase (light blue squares) and the fully polarized normal 
(FPN) phase (dark blue diamonds). The presence of the ID optical potential leads to a significant 
enlargement of the FFLO region. We used s = ^j - = 5 and A = 1200nm for the depth and wave 
length of the optical potential, and \a s \ = 500nm for the s-wave scattering length. 



is determined for a discrete set of (/i, C)-points and each of these points is displayed. We 
chose this type of visualization so that later on in this paper, the transformation to the phase 
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diagram at fixed densities will be clarified by the transformation of the set of (//, C)-points 
to a set of Cj^-, P^j -points, where P denotes the polarization. In both phase diagrams in 
figure [TJ the system is always in a spin-balanced superfluid state when ( = 0. This is not 
visible in Fig. [1] since in this figure ( > 0.06. For small values of the total chemical potential 
fi (fi < 0.5 in Figs. ED(A) and (B)), the system makes a transition into a partially polarized 
normal state when £ becomes of the order of the gap A. For larger values of the system 
goes through all three different states when the imbalance chemical potential ( is increased 
from zero. For low values of the system remains a spin-balanced superfluid. This means 
that a finite energy is needed in order to convert particles from one spin species to the other. 
The minimal amount of energy that is required for this conversion increases when the total 
chemical potential increases. For high enough values of £ the Fermi gas becomes polarized. 
At this point the system makes a transition into the FFLO state, both in figure [I] (A) and 
(B). Characteristic of the FFLO state is that it has an oscillating order parameter 23j, \24 \. 
This order parameter is enhanced by the presence of the ID potential, which results in a 
much larger FFLO region in Fig. [1] (B) when compared to the case without the optical 
potential in Fig. [H(A). When the imbalance chemical potential ( increases further, the sys- 
tem will go over into a partially polarized (PPN) state and eventually into a fully polarized 
(FPN) state, and superfluidity is destroyed. 

In the phase diagram in Fig. (TJ we made the following choice for the depth s = Vq/Er 
and for the wave length A of the optical potential: s = 5 and A = lOOOnm, which are 
both realizable in current experiments. The depth of the optical potential must be chosen 
carefully, because the dispersion relation ([3]) that is used in our description is only valid in 
the tight-binding limit. This will be discussed in more detail in Sec. IIII A[ Both parameters 
s and A correspond to specific values of 6 and Ql, which are defined in expression ([3]). Since 
we use units of scattering length, a fixed value of this quantity has to be chosen, in order 
to compute the values of S and Ql- Here we choose \a s \ = 500nm, which corresponds to a 
situation close to a Feshbach resonance. 



B. Phase diagram at fixed densities 

In the previous section, we chose to construct a phase diagram at fixed chemical potentials 
fi and (. This is not the only possible choice. Another possibility is to work at fixed total 



7 



density n = + nj, and fixed polarization P = (n-j- — n_\) / (n-j- + nj.), which is equivalent 
to fixing the density of both spin species and n±. To transform between these two 
descriptions, the following number equations have to be solved: 

dQ sp (/i, C; A, Q) 



n 



5n 



dQ sp O, C; A, Q) 



(5) 



In principle, we need to use the chain rule in since A and Q are separate variables that 
depend on /i and (. But the saddle-point condition (J4j) will make these terms disappear. If 
fluctuations around the saddle point are taken into account, these extra terms will remain. 
The inclusion of fluctuations around the saddle point lies beyond the scope of this paper. 
Starting from the phase diagram at fixed chemical potentials in figured] (B), the phase dia- 
gram at fixed densities can now be constructed. Every point in the phase diagram depicted 
in figured] (B) corresponds to a set of values for the parameters (fi, £; A, Q). Inserting these 
values into both equations in fl5]) yields the density n and the density difference 5n = n^—n±. 
The latter can also be written as the polarization P = 5n/n. The final step is to find the in- 
teraction strength (/cpa s ) -1 that corresponds to the density n. Since the scattering length a s 
is kept constant, this can be done by using the following expression that relates the density 
to the Fermi wave vector k F - 
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This expression ([6]) is derived by calculating the number of particles at temperature zero 
that have an energy that lies lower than the Fermi energy, where the particles are assumed 
to have a dispersion relation given by (]3]). The procedure outlined above yields the phase 
diagram at fixed density and at fixed polarization, which is shown in figure [2] Each point 
on this phase diagram corresponds to a point on the (/z, C)-pliase diagram in figured] (B). In 
the case of zero polarization, the system is in the spin-balanced superfluid state. Figure [2] 
shows that there exists a portion of the phase diagram at finite polarization, where no points 
are present. Mathematically, this corresponds to the case when there exists no combination 
of density and polarization that can satisfy both number equations ([5]). Physically, this is 
a state where the system separates into two different phases: a spin-balanced superfluid 
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(k F a s )"' 

FIG. 2: Phase diagram of an imbalanced Fermi gas in 3D, subjected to a ID optical potential, 
at fixed density and fixed polarization P. The dimensionless parameter (fc^a s ) _1 characterizes the 
interaction strength. Each point in this phase diagram is a transformation of a corresponding 
point in the phase diagram in figure [I] (B). At polarization zero, the system is in the spin-balanced 
(S-b) superfluid phase. When polarization becomes finite, there is a region where no points exist. 
This means that no combination of density and polarization can solve the number equations, 
and the system separates into two phases: a spin-balanced superfluid component and an FFLO 
component. Interestingly, both the boundary between the phase separation region and FFLO, 
and between FFLO and the normal phase have a maximum (indicated by a red square and a red 
circle respectively) . The dashed oblique line indicates the boundary where the Fermi energy of the 
majority spins becomes equal to the top of the first Bloch band, which is equal to 25. The empty 
blue squares (normal phase) and empty green triangles (FFLO) indicate that the Fermi energy 
of the majority spins lies higher than the top of the first Bloch band. The vertical dashed line 
indicates the value of the interaction strength where the Fermi energy at polarization zero becomes 
equal to the top of the first Bloch band. Both dashed lines intersect with one of the two indicated 
maxima, (s = 5, A = 1200nm, and \a s \ = 500nm) 

phase and an FFLO phase, where the latter is being pushed out of the former. When 
polarization increases further, the system makes a transition into the FFLO state. Figure 
|2] also shows that when the interaction strength increases, both the boundary between the 
phase separation region and the FFLO state and the boundary between the FFLO state 
and the normal state, lie at higher values of polarization. Interestingly though, for both 
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boundaries there exists a critical interaction strength at which this polarization becomes 
maximal. These two maxima are indicated by a red square and a red circle respectively in 
figure If the interaction strength becomes larger than these two critical values, both the 
region of phase separation and the FFLO region seem to shrink. This counterintuitive effect 
will be discussed in the next section. When polarization increases further, the energy cost 
of forming pairs with finite momentum will eventually become too high and the system will 
make a transition into a normal Fermi gas. The points with polarization P = 1 correspond 
to the fully polarized region in Fig. [TJ We now have given a general overview of the phase 
diagrams of an imbalanced 3D Fermi gas subjected to a ID potential. In the next section 
we will investigate the effects of the properties of the optical potential on the FFLO state. 

III. EFFECTS OF THE ID POTENTIAL ON THE FFLO STATE 

A. Role of the wave length and the depth of the optical potential 

There are two parameters that determine the properties of the ID optical potential: the 
potential depth s, and the wave length A. Before discussing the effect of these parameters, 
it is important to know to which intervals they are limited in our current treatment. For 
the wave length, we take realistic values corresponding to present day experiments: A G 
[600nm, 1200nm]. We have to be careful, however, when choosing values for the depth of 
the optical potential. There are two physical reasons to set a lower boundary to this depth. 
Firstly, since the optical potential is described by a tight-binding dispersion relation (Eq. [3]), 
the depth of the potential must be large enough to justify this approximation. To quantify 
this, we calculated the energy spectrum of a particle in a periodic potential, and compared it 
to the tight-binding dispersion given by ((3]). For completeness, we added a concise overview 
of this calculation in appendix [A] This calculation shows that the error made by using 
expression ([3]) for the dispersion is less than five percent when s > 4. Secondly, by using the 
dispersion relation ([3]), we only take into account the lowest Bloch band. This description 
is justified as long as the system has a Fermi energy which lies below the second Bloch 
band. A Fermi energy that lies between the top of the first Bloch band and the bottom of 
the second Bloch band is allowed, since the dispersion in the direction perpendicular to the 
direction of the optical potential is still the same as in the free 3D Fermi gas case. We have 
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verified that by considering a potential depth s > 4, the system has a Fermi energy that lies 
below the bottom of the second Bloch band, for all cases considered in this paper. Based 
on this analysis, we will use the lower boundary s = 4 for the potential depth throughout 
the rest of this paper. There is also a reason to set an upper limit to the depth of the 
optical potential. When the potential becomes too deep, the 3D Fermi gas will transition 
into a series of weakly-coupled two-dimensional (2D) pancakes. Such a 2D Fermi gas must 



be described using a different renormalized interaction strength [35(. We therefore choose s 
to lie between 4 and 6, so that both the 3D description and the tight-binding approximation 
are justified. 

We now proceed to investigate the effects of both optical potential parameters A and s. 
Figure [3] shows several phase diagrams at fixed density and polarization, for an imbalanced 
Fermi gas in 3D subjected to a ID potential, for various values of both parameters A and s. 
In Figs. [3] (A), (B), and (C) the wave length of the optical potential is varied. These three 
figures show that the maximal polarization at which the FFLO state can exist increases 
when the wave length of the optical potential becomes smaller. This effect agrees with our 
expectations and can best be explained in terms of the Fermi surfaces in momentum space. 
When an imbalance exists between spin-up and spin-down particles, an energy gap opens up 
between the Fermi surfaces of the two spin species. This is akin to a Zeeman-splitting effect. 
In order to bridge this gap and to let particles pair with the same energy, the FFLO state 
forms pairs with a finite momentum Q. The larger the imbalance, the larger the energy gap 
becomes and the larger Q has to be. In Ref. 30] we demonstrated that the wave vector of 
the FFLO state cannot become larger than the wave vector of the optical potential. Since 
a smaller wave length corresponds to a larger wave vector, the maximal value of the wave 
vector Q that can be attained by the FFLO state increases when the wave length of the 
optical potential decreases. This means that a larger energy gap can be bridged and thus a 
larger polarization can be accommodated. Figures [3] (C) and [3] (D) show the phase diagram 
of the system at the same wave length but at different potential depth. In both figures, the 
maximal polarization of the FFLO state is equal. We conclude that only the wave length 
of the optical potential has an influence on the capability of the system to remain in the 
FFLO state while being polarized, at least within the small range of potential depth s that 
is considered here. 
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FIG. 3: Overview of several phase diagrams of an imbalanced Fermi gas in 3D subjected to a ID 
optical potential, for different values of the depth s and the wave length A of the potential. The 
maximal polarization P at which the FFLO state can exist (indicated by a dashed horizontal line) 
increases when the wavelength of the optical potential becomes smaller. At constant wave length, 
the depth of the optical potential does not have an influence on this maximal polarization, at least 
within the range 4 < s < 6 that is considered here. (S-b SF = spin-balanced superfluid, PS = 
phase separation) 

B. Effect of the edge of the Brillouin zone 

One effect that is visible throughout all phase diagrams in figures [2] and [31 is that when 
the interaction strength exceeds a certain critical value, both the region of phase separation 
and the FFLO region become less stable with respect to polarization (i.e. the maximal 
polarization at which they can exist, decreases) . This can be seen by the fact that both 
the boundary between the phase separation region and the FFLO state and the boundary 
between the FFLO state and the normal state have a maximum (indicated by a red square 
and a red circle respectively in figure [2]). Moreover, these two maxima seem to correspond 
to different values of the interaction strength. This effect is counterintuitive, because when 
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the interaction strength increases, the binding energy of the fermion pairs increases likewise, 
which in general should strengthen superfluidity, not weaken it. 

This effect is explained with the help of figure HJ In this figure, the empty blue squares 
(normal state) and the empty green triangles (FFLO state) indicate that the Fermi energy 
of the majority spins lies higher than the top of the first Bloch band, which is equal to 28. 
The boundary at which the Fermi energy of the majority spins becomes equal to the top 
of the first Bloch band is indicated by the dashed oblique line. The maximum polarization 
under which the FFLO state can exist (red circle), lies exactly at this line. Furthermore, the 
vertical dashed line indicates the value of the interaction strength at which the Fermi energy 
at polarization zero becomes equal to the top of the first Bloch band. This dashed vertical 
line intersects with the maximum polarization under which the phase separation region can 
exist (red square). The reason why each of these two dashed lines intersects with one of the 
two maxima can be explained as follows. Consider the case which is shown schematically 
in Fig. HJ Figure H] (A) shows the Fermi surface of both spin-up and spin-down particles 
at polarization zero. In order for the system to become superfluid, pairing occurs between 
particles that lie in a band with thickness A around the Fermi energy. When these pairing 
bands of the spin-up and spin-down particles overlap, pairing is possible and the system is 
in the superfluid state. When polarization is introduced, a gap arises between the two Fermi 
surfaces, as indicated in figure|4](B). As long as this gap is smaller than A, there is overlap of 
the pairing bands, and superfluidity is possible in principle. When the polarization becomes 
too large, the overlap disappears and BCS superfluidity is no longer possible. In order to 
remain superfluid, the system can transition into the FFLO state, which translates the Fermi 
surface of the minority spins so that it locally lies on the Fermi surface of the majority spins, 
making pairing possible again (see Fig. H](B)). This allows the system to remain superfluid, 
at the cost of paying some energy to form particles with finite center-of-mass momentum. 
This is the case for a 3D Fermi gas, but when a ID potential is added, an additional effect 
comes into play. When the density and/or polarization is high enough so that the Fermi 
surface of the majority spins lies at the edge of the first Brillouin zone (this is true for the 
area at the right side of the dashed oblique line in Fig. [2]), pairing at this Fermi surface 
becomes more difficult, because for k z > Ql, there are no available states (see Fig. @] (C)). 
Hence, the formation of an FFLO-type superfluid phase will be hindered. This now explains 
why, at equal polarization, the normal phase becomes favored over the FFLO phase when 
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the interaction strength exceeds the value corresponding to the situation where the Fermi 
energy of the majority spins equals the top of the first Bloch band (indicated by a red circle 
in Fig. [2]). The question can be raised whether the FFLO state can still exist if the fermion 
pairs acquire a momentum in the x or y direction, since there exists no Brillouin zone in 
these directions. We did not treat this case in our paper because we assumed the FFLO 
momentum to be of the form Q = (0, 0, Q z ). 

In the case of figure H] (C), BCS superfluidity can still exist, because when the Fermi 
surface of the minority spins is not translated, pairing does not occur at the boundary of 
the first Brillouin zone. However, when the density becomes high enough, even the Fermi 
surface of the particles at zero polarization will lie at this edge. This situation is depicted in 
figure H] (D). For this and for higher densities, pairing of BCS superfluidity will also occur 
at the edge of the Brillouin zone and it will suffer the same problem as the FFLO state. 
This effect can be seen in Fig. [2] from the fact that, at equal polarization, the FFLO state 
becomes favored over the spin-balanced superfluid phase, when the value of the interaction 
strength increases above the value corresponding to the situation where the Fermi energy 
at zero polarization equals the top of the first Bloch band (indicated by a red square in Fig. 
E]). 

IV. CONCLUSIONS 

We have investigated the effects of a ID optical potential on the FFLO state in an 
imbalanced Fermi gas in 3D. To allow for a direct connection with the experiments, we 
constructed the phase diagram of the system, both at fixed chemical potentials and at fixed 
densities of both spin species. Subsequently, we investigated the effects of the depth and the 
wave length of the optical potential on the FFLO state. It was shown that this state can exist 
under a larger polarization, when the wave length of the optical potential becomes smaller. 
The depth of the optical potential did not have any appreciable effect, within the range 
that can be considered in our tight-binding approximation. By studying the phase diagram 
at fixed density we have discovered an unexpected effect. When the interaction strength 
exceeds a certain critical value, the maximal polarization at which both the FFLO phase 
and the phase separation region can exist decreases. The underlying reason was discussed in 
detail. We have found that when the density of the system becomes so high that the Fermi 
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FIG. 4: Schematic representation of the pairing mechanisms in an imbalanced 3D Fermi gas sub- 
jected to a ID optical potential. For the sake of clarity, we have omitted the y direction. (A) The 
Fermi surfaces of the spin-up and the spin-down particles at polarization zero. Pairing occurs in 
a band with thickness A around the Fermi energy and BCS superfluidity is possible. (B) When 
polarization is introduced, a gap is created between the Fermi surfaces of both spin species. If this 
gap becomes too big, there will be no overlap between the pairings bands around the Fermi sur- 
faces. The system can solve this by translating the minority Fermi surface (dashed red ellipsoid) 
so that locally, both Fermi surfaces are re-aligned (indicated by a green rectangle). (C) When 
the density increases, the Fermi surface of the majority spins (here spin-up) will lie at the edge 
of the first Brillouin zone. This will hinder FFLO-like pairing because no states are available at 
k z > Ql- (D) Above a critical density, the Fermi surface lies at the edge of the Brillouin zone even 
for polarization zero. At this point, BCS superfluidity will also be hindered. 



energy exceeds the top of the lowest Bloch band, superfluid pairing becomes hindered. This 
is because in that case, particles with the highest energy lie at the edge of the first Brillouin 
zone. Since no states exist outside of this zone, pairing is not possible and the superfluid 
state is suppressed. We expect this effect, along with the enhancement of the FFLO state, 
to be observable in experiments where a ID optical potential is applied to an imbalanced 



15 



Fermi gas in 3D. 
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Appendix A: Range of optical potential depths that allow the tight-binding approx- 
imation 

The energy spectrum of a particle in a periodic potential can be calculated using 
Schrodinger's equation. The potential under consideration is given by 

V(z) = sE R sm 2 (Q L z), (Al) 

where s is the depth of the optical potential in units of the recoil energy E R and Ql is the 
wave vector of the optical potential given by Ql = ir/\. Given this potential, Schrodinger's 
equation states that 

h 2 d 2 fl-cos(2Q L z) 



2m dz 2 R 



i>k (z) = e k ip k (z) . (A2) 



Since we consider a particle in a periodic potential, Bloch's theorem applies and we can 
re-write the wave function ip k (z) as 



*p k (z) = e lkz u k (z) , (A3) 

where u k (z) has the same periodicity as the potential. Because u k (z) is periodic, it can be 
written in a Fourier series 

1>k {z) = e^ (X>e^ 

= J2cne l(k+2nQL)z . (A4) 
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Substituting (|A4]) in Eq. fl£2]) yields 
r ^2 



E 



i{k+2nQ L )z 



sE f 



E 



D 2iQ L z 



+ e 



-2iQ L z 



i(k+2nQ L )z 



E' 



£fc , c n e 



i{k+2nQ L )z 



(A5) 



Re-indexing the second sum over n and identifying the coefficients on both sides in equation 
i5l) leads to 



' h 2 

— (k + 2nQ L ) 2 
2m 



sE 



R 



sE 



R 



Cn+1 + Cn-l, 



0. 



(A6) 



This expression represents an infinite set of linear equations. In practice, it is enough to 
consider five equations, since the contributions of the eigenvectors c n decreases exponentially 
for large n. Solving equation flA6j) numerically yields the exact energy spectrum of a particle 
in a periodic potential, given by (1A1I) . Figure [5] shows a comparison between this exact 
spectrum and the corresponding dispersion relation in the tight-binding limit, given by Q. 
For s < 2 the difference is substantial, but for s > 4, the approximation deviates less than 5 
percent from the exact result. This justifies the approximation that we make by considering 
s > 4. 
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